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Hidden Symmetries in Nano-magnets 
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The hidden symmetry of certain nano-magnets leads to many of the levels being doubly degenerate 
for periodic values of the Zeeman energy. Corresponding to such a symmetry is an operator, K n , 
related to the time reversal operator, and which commutes with the Hamiltonian. The degeneracies 
for whole integer spin values are often similar to the Kramer's degeneracy for half-integer spin. 
An all algebraic method for constructing Hamiltonians with such hidden symmetries and different 
crystal symmetries is described. 

PACS numbers: 75.45+j,75.50Tt,75.60Ej 



For both of the much studied molecular magnets 
Mni2H and FegM it is found experimentally that most 
of the levels cross, thereby becoming doubly degenerate, 
at a set of periodic values of the magnetic fields applied 
along a certain symmetry direction. In fact, neither of 
these statements is precisely true. Levels cross at almost 
the same field and the fields at which this occurs are 
almost periodic up to some maximum value. Theoret- 
ically this phenomena of simultaneous degeneracies has 
been described in terms of topological interference at "di- 
abolic" points|| and in terms of intermediate spinQ. In 
both these theoretical approaches certain approximations 
need to be made. 

In this Letter it is shown that the degeneracy of nearly 
all the levels reflects a hidden symmetry, i.e., that there 
exists an operator, K n , related to the time reversal oper- 
ator which commutes with the Hamiltonian for the spe- 
cial values of the field B n . It will thereby be shown how 
to construct spin models which have the property that 
many of the levels cross at exactly the same field for ex- 
actly periodic fields (up to some maximum), this inde- 
pendent of the spin value S and even if the parameters 
correspond to the tunneling regime. In effect, it is shown 
that whole integer spins, in a suitable finite magnetic 
field Bi, exhibit a Kramer's degeneracy, for all but a 
single highest energy level. This supports the intermedi- 
ate spin interpretation^] of the situation since in effect a 
whole integer spin has been converted into one which for 
practical purposes behaves as if it was half integer. Un- 
like other approaches to this problem the present one is 
purely algebraic and without approximations. (Although 
the same algebraic methods will be used to generate use- 
ful approximate results.) The result has implications for 
magnetization plateaus for concentrated bulk magnets. 

Experimentally relaxation via tunneling is exactly sup- 
pressed at applied fields for which there is such a hidden 
symmetry. This has implications for applications in both 
magnetic recording and quantum computing. 

The large spin model appropriate to Feg is traditionally 
written p| as, 



where without loss of generality the positive parameters 
have D > and where the coordinate system coincides 
with the principal crystal axes. 

The first exercise is to show for this model all the lev- 
els but one exhibit a Kramer's degeneracy when B = 
h x j2E(D + E). The existence of this massive degener- 
acy is intimately related to the structure of the SU(2) Lie 
algebra. Consider the most general linear combination 
of the elements of this algebra, written as, 



r = aS x + ipSy, 



(2) 



where S z has been eliminated by a suitable choice of axes. 
Take T to be the time reversal operator for integer spin, 
so T 2 = 1, then, 



-TrT, 



(3) 



independent of the c-numbers a and (3. The pseudo-time- 
reversal-operator is, written as, 

Ki = Tr. (4) 

This has the evident property that [H, K{\ = 0, if, 



Using the definition of r with a = 
V2E, 



(5) 



x/D + E and /3 



Hi = -(£> + E)S X 2 - 2ES y 2 + hy/2E{D + E)S Z , (6) 

which, to within a constant, is equivalent to Eqn. (jl]) 
with gusB = hyj2E(D + E) and with a change of axis 
so that x is the easy direction. If \n) is an eigenstate 
of Hi with energy E n then K\\n) is either (i) zero or 
(ii) another degenerate eigenstate of Hi. In the stan- 
dard way, HiK\\n) = KiHi\n) = E n Ki\n), then for the 
normalization, 



(n| r^T^Tr \n) = (n\ rV |n) 



-En 



(7) 



H 



-DSy 



E[S X 



gH B BS x 



(1) 



since T is Hermitian, i.e., T^T = 1 and using Hi = —r^r. 
Thus (i) if K\ destroys \n) then E n — and otherwise 
(ii) when the norm of Ki\n) is finite then E n < 0. It is 
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necessary to show that Ki\n) and 
To this end it is useful to write: 



are different states. 



(b) 



Wi 



-r'r 



-(a 2 - (3 2 )(S+ 2 + S- 2 ) 
-(a + (3) 2 S-S+ - (a- [3) 2 S + S~, (8) 



which makes evident the fact that Tt breaks into two 
blocks. It is possible to insist that there exists a nor- 
malized \n, e) = J2m fl 2m I S, 2m) which involves only 
the eigenstates of S z with even eigenvalues 2m, i.c, 
S z \S, 2m) = 2mh \ S, 2m). Similarly there is the normal- 
ized odd state |n, o) = a2m+i \S, 2m +1). It is then 
observed that K x =Tr = T[(a + /3)S~ + (a - f3)S+], so 
that K\ \n,e) oc \n,o) and \n, e) and \n,o) are indepen- 
dent. 

Since for whole integer S, the total number of states, 
2S + 1 , is odd there must be at least one state which is 
not part of a doublet and therefore for which E n = 0. 
That this state is unique is proved by explicit solution of 
the equation K\ |n) = 0. The resulting matrix equation 
is, 
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= 0. 



(9) 



\ a-s+i sinaM~f +1 / 

where M™ = [S(S + 1) - nm] 1 ' 2 . This reduces to: 



os-i = a S-3 
«S-2 : 



sinaAff 1 
5^2 a s 



cos aM 



a-s+i = 
as-i = 



and 
sin aM 



s-i 
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sin aM-s+2 
cos aM 



-a-s+2- 



S-3 
5 — 2 

snas-2, 
(10) 



-s+i 



The solution is uniquely determined up to a multiplica- 
tive constant, proving the uniqueness of the state with 
E n =0. 

The situation is illustrated in Fig. |l|a, all the lowest lev- 
els form pseudo-Kramer's doublets with the sole excep- 
tion of the highest energy state which has zero energy. 
If one ignores the highest energy singlet, the spectrum 
is identical to that for half- integer spin S = S — (1/2). 
These conclusions are independent of the parameters D 
and E and the value of S, i.e., these precise degeneracies 
occur even in the absence of a tunneling barrier. 



FIG. 1: (a) Illustrated for S — 4 is the situation for the first 
magic field gusBi — Yi\j2E(D + E). There are four pseudo- 
Kramer's doublets with a highest energy singlet. The singlet 
is at zero energy, (b) At the second magic field B2 for S = 4, 
there are three doublets with three singlets at higher energy. 



It is the case that [r, r f ] = 2hyj2E(D + E)S Z and so 
the second Hamiltonian in the series is, 



H2 = —r^r + [r, r 1 '], 



which corresponds to g[ i B B 2 = 3hy/2E(D + E). This 
also has pseudo-Kramer's degeneracies. To prove this it 
is necessary to find a K2 such that [H2, K2] — 0. Since 
r is related to the raising operator a first guess might be 
Tr 3 , however the relevant commutator is rather, 

[H 2 , Tr 3 ] = 2af3(/3 2 - a 2 ) [Tr, S z ] ? 0. (11) 

The commutator to the right is that involved in 
[W 2) ifi] = -2a{3[Tr,S z ], and it follows that, 



K 2 = Tr 3 + ((3 2 - a 2 )Tr. 



(12) 



The special case E = D, when H 2 = 2ES z (S z + 3) and 
K2 = T(S+) 3 , then establishes that K 2 destroys three 
states which are sing lets, Fig. |b. 

The generalization is straightforward: TL n = —r*r + 
(n — l)[r, r'] and K n has the form: 



Kn 



2\n-k-l 



Tr 



2k+l 



(13) 



fe=0 



where a\ are integers, and where it is possible to set 
a" -1 = 1. The power of (ji 2 — a 2 ) ensures that all terms 
(f3 2 — a 2 ) n - k - 1 Tr 2k+1 have the dimension (Energy)"" 1 " 2. 

It is an interesting but tedious exercise in linear alge- 
bra to construct an algorithm to determine the a\. The 
results for small n are: 



n = 3 
n — 4 



al = 5, a§ = 4, 



14, 



49, al = 36. (14) 



The method and results for several other values of n are 
given elsewhere ||. Again, the special case when E = D 
illustrates that there are 2n — 1 singlets for a field B n = 
(2n — \)h^/2E(D + E) so that n = S corresponds to the 
largest value of B n for which there is a pseudo-Kramer's 
doublet. 

In summary, it has been proved that there exist a series 
of operators K n which commute with the Ji of Eqn. (|l|) 
for exactly periodic magic values B n of the applied field 
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along hardest axis of the nano-magnet. For such values of 
the field there are 25 + 1 — 2n pseudo-Kramer's doublets 
which are exactly degenerate. 

The highest symmetry of Eqn. (|l|) is two-fold. It is 
interesting to investigate the problem with a four-fold 
axis. The simplest possibility is to set, 

r = f]S +2 +jS- 2 (15) 

whence r* — TrT. As before define K — Tr, then Hi = 
—K 2 can be written as, 

Hi = -K 2 = -( v 2 + 7 2 )[5 Z 4 + (5 - 25(5 + l))S z 2 } 

-{r, 2 - 7 2 )[45 z 3 - (2 - 45(5 + l))S z }. (16) 

This H has the desired symmetry. The diagonal 
anisotropy term has the ratio of the coefficients of S z 
and S 2 fixed and the Zeeman term contains S z 3 in ad- 
dition to S z . For all values of 5 this corresponds to an 
easy plane magnet. Notice that the off-diagonal part 
connects states with S z quantum numbers 

which differ by four and implies that the matrix H breaks 
down into four blocks. Consider whole integer spin 5, the 
blocks behave differently as the values of S z are even or 
odd. The operator K acting upon a solution, \n) which 
contains the state with S z = produces a state K\n) in a 
different block, i.e., that which involves the S z = 2 state. 
However when K operates on a state which contains ei- 
ther of the S z = ±1 states, the result belongs to the same 
block and the limit 777 — > shows K\n) is not distinct 
from |n). This should be compared with the situation 
when the Zeeman term, (r] 2 -j 2 )[AS z 3 ~(2-AS(S+l))S z }, 
is absent. Then the plain time reversal operator T com- 
mutes with H and connects the S z = +1 to the S z = — 1 
blocks but now maps the solutions in the S z = and 
S z = 2 blocks back to themselves. Thus, e.g., if 25+ 1 is 
a multiple of four there are half of the states which are 
degenerate both in the absence of a Zeeman splitting and 
for the sequence of Hamiltonians H n = Hi + (n — 1) [r, r^] 
and operators K n constructed by the present methods, 
i.e., these H n correspond to the degeneracies of a whole 
integer point rather that a half integer point for which 
Kramer's theorem insists that all states form doublets. 
A second possibility is to take, 

r = ^ 1/2 (5+ + -^5- 3 ), (17) 

whence = —TrT and, 

Hi = K 2 = (Tr) 2 = -rV 

= -DS-S+ - a (5+ 4 + 5- 4 ) - ^5+V 3 

= DS 2 - a (V 4 + 5~ 4 ) 

+H+(D + %-Z)S z , (18) 



where H = -DS{S+l)-{a 2 /D)[M 4 -8M 2 + M(23S z 2 + 
12) - 285 z 2 , with M = 5(5 + 1) - S 2 , and where the 
Zeeman term involves, 

Z{S 2 )= 3 [25 2 (5 + l) 2 + 125(5+1) + 4] 

- [185(5+1) - 51]5 Z 2 + 95 z 4 . (19) 

With this choice of r is is possible to change the ratio of 
5 2 4 and S z 2 in the diagonal anisotropy but at the cost of 
a more complicated Zeeman term, [r, r'] — (D+ jyZ)S z . 
However both terms simplify in the limit that D a 
whence the magic values have g^sBn = (2n — l)hD. 
The operators K n now connect even and odd S z blocks 
and so the H n = Hi + (n — l)[r, r>] and K n correspond 
to the double degeneracies of a half integer point as for 
the development with two fold symmetry. 

As an alternative approach to four fold symmetry, con- 
sider the Hi, of Eqn. (§) with E = 1 and D = 3. This 
can be written as, 

Hi = -r^r = -25(5 + 1) + 2(5 Z 2 - S x 2 ) + 2V2S Z . 

(20) 

The square of Hi + 25(5 + 1) - 1 = -1 + 2{S 2 -S x 2 ) + 
2V2S Z is, 

7U = [fti + 25(5+l)-l] 2 

= l + 4(5 z 2 -5, 2 ) 2 + 4[5(5+l)-5 ?/ 2 ] 

+8x/2(5 z 2 - S X 2 )S Z , (21) 

which in the absence of the non-trivial Zeeman term 
8\/2(S z 2 — S X 2 )S Z has x «-> z symmetry. In generaliz- 
ing this approach to n =/= 1 this Zeeman term will change 
form with n. 

The conclusion reached by all three methods is that 
the occurrence of simultaneous degeneracies with a four 
fold symmetry implies a non-trivial Zeeman term. In 
nature molecular magnets can and will have such non- 
trivial time reversal symmetry breaking terms. For ex- 
ample, when there is finite angular momentum in the 
crystal field ground state, and when this corresponds to 
a multiplet which is reflected by a fictitious spin, then 
it is entirely possible to have complicated Zeeman terms 
of the type which arise here. Probably the best known 
example is the Tg ground state of rare earth ions in cubic 
symmetry. 

Although the simplest model with four fold symmetry, 

H = DS 2 - a (s +i + 5~ 4 ) + gfi B BS z , does not have 

exact pseudo-Kramer's degeneracies, the present formal- 
ism can be used to determine approximately the fields at 
which the lowest four levels cross for not too small 5. 
For positive D > a, the models, Eqn. ( |l6| ) and Eqn. (p^), 
yield different information. With this sign of D the diag- 
onal anisotropy term is dominated by the S z 2 term and 
this localizes the wave function near S z = 0. For model 
Eqn. (|B|), 25 2 (?7 2 + 7 2 ) = D and the Zeeman term is 
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in turn dominated by AS 2 (r] 2 — r y 2 )S z « 2DS Z . Directly, 



for model Eqn. (18) the coefficient of S z is D while the 



dominant part of the Zeeman term is only DS Z . The 
difference in the two Zeeman terms reflects the different 
nature of the degeneracies as explained above. For inte- 
ger spin S, the field such that g/isB = HD deduced from 
Eqn. (|l^) corresponds to the first half-integer point with 
two doublets while g^BB = 2hD deduced from Eqn. ( pjj ) 
corresponds to the first finite field whole-integer point 
with a doublet and two singlets (and a singlet ground 
state). 

The model Eqn. (|l^) can also be used to estimate devi- 
ations from periodicity. Again consider specifically posi- 
tive D > a and modestly large whole integer S. The low- 
est lying states localize near S z ~ g/isB /2hD which im- 
plies an approximation: Z(S Z 2 ) => Z((g/iBB /2hhD) 2 ). 
Then e.g., for the low lying levels the first half-integer 
field B\ is approximately determined by solving, 



g t x B B 1 = (D + —Z((g^ B B 1 /2hD) 2 )). 



(22) 



It might also be noted that the anti-linear nature of 
T is not essential to the development. Consider some 
other "reflection" symmetry R such that R 2 = 1. This 
must conserve the angular moment commutation rules 
so some components of S change sign while others do 
not. For example, there exists a R such that RS X R = 
—S x , RS Z R — — S z but RS y R = S y for which RrR = 
—r^ and for which the above development is essentially 
unchanged. 

In the much used functional integral method degen- 
eracies correspond to the destructive interference of two 
equivalent tunneling paths and this implies that the topo- 
logical phases are ±7r/2. This can also be verified us- 
ing the current analytic techniques. Consider again the 
model Eqn. @, the topological phase can be "mea- 
sured" as a Berry phase || for a suitably conceived adia- 
batic process. The two states \n, a), a — e or o have 
(S x ) = 0. The imagined adiabatic process is there- 
fore between the states |n, ±) = (2)~ 1 / 2 [|n, e) ± |n, o)] 
which have (S x ) ~ ±5*. The tunneling path which passes 
via the positive or negative y-axis is favored by adding 
^fhtSy, h t = gj-isBy > 0. Observe first that with h t — 
the matrix TL n is real and therefore solutions |n, ±) can 
be made real. On the other hand S y is a pure imaginary 
matrix. It is trivial that there are pure imaginary ma- 
trix elements proportional to h t which connect the |n, ±). 
With a very small h t the magnetization will therefore 
slowly and approximately adiabatically oscillate between 
S x ~ ±S with one of the above least action paths domi- 
nating. Reflecting the pure imaginary matrix elements, if 



the initial state is \n, +) then after a half-period the sys- 
tem will find itself in the state — zsign(/i t ) \n, —) demon- 
strating that the Berry phase is = F 7r /2- 

Especially in one dimension, the nature of the ground 
state of a concentrated magnet depends on the half or 
whole-integer nature of the spins with whole integers 
having gaps|Q. It is possible to imagine, e.g., a mag- 
net with Eqn. (Q) describing the intra-atomic interactions 
and with, e.g., — Y)^ JijSi-Sj interatomic interactions. If 
the spin S is whole integer, and if J is suitably small com- 
pared with the intra-atomic splittings ~ D, then a field 
with one of the magic values B n will transmute the sys- 
tem to one with effective half-integer spin S— (2n— l)/2. 
Since the B n values are periodic the phase diagram will 
also be periodic in the applied field. 

The period measured in terms of the magnetization per 
site can be determined assuming that the magnetization 
is given to a good approximation by the classical result. 
(Even for a small spin value, say S — 2, this is a good 
approximation for most parameters ranges.) The term 
— (D + E)S X 2 might be considered to be a x-directed 
anisotropy field Ba such that g/isBA = 2Sh(D + E). 
In the presence of a real field B along the z-direction 
the dimensionless magnetization component along this 
direction is m w S t&n{h/2Sh(D + E)) « (h/2h(D + E)) 
and the period of the whole integer points is Am = 
\J Ef (D + E) < I. Thus magnetization plateaus associ- 
ated with the gaped whole integer phases will not occur 
at simple multiples of gfiB as is often supposed to be the 
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